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Abstract. In this article we consider the Boussinesq system supple- 
mented with some dissipation terms. These equations model the prop- 
agation of a waterwave in shallow water. We prove the existence of a 
global smooth attractor for the corresponding dynamical system. 

1. Introduction 

This article is concerned with the long time behavior of the solutions to a 
damped-forced Boussinesq system that read 

f nt + u x + (r}u) x - rj xx = 0, 

\ Uf - U txx ~ Uxx + Vx + u x u = f. 

Here we have an incompressible fluid on a channel. u(t, x) is the horizontal 
velocity at the top of the fluid, rj is the fluctuation of the height of the fluid 
with respect to the rest position that is z = rj(t, x) = 0, assuming that the 
bottom of the channel is at z = — 1. Observe that in our model we have to 
ensure that r](t,x) > —1 Vi,x. 

Here f(x) is an external force that does not depend on time and the 
damping terms are respectively —u xx ,—rj xx . In the conservative case, that 
read 

f Vt + u x + (r]u)x = 0, q\ 

\ u t - utxx + Vx + u x u = 0; 

this system has been introduced by Boussinesq in 1877 to model the fluctua- 
tion of a waterwave in shallow water. Other well-known asymptotic models 
are Korteweg-de Vries equations and Benjamin-Bona-Mahony equation, also 
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known as the regularized long wave equation. For these asymptotical models 
we would like to refer to [5, 19] and to the references therein. 

In this article we are interested in the dissipative case. In the case where 
/ = 0, the solutions converge to the equilibrium and the issue is to find out 
the rate of convergence. Following the pioneering work of Amick, Bona and 
Schonbek [4], this issue has been addressed in the case where x G M. D , D > 2 
using the famous Schonbek splitting method [16]. Here we plan to study the 
dynamical system provided by (J2J into the framework of infinite dimensional 
dynamical system [11,15,17]. Our main result states as follows 
Theorem 1.1. The dynamical system provided by £JJ^ features a compact 
global attractor into a suitable energy space. Moreover this compact global 
attractor has finite fractal and Hausdorff dimension. 

This result compares with previous results obtained for dissipative KdV 
equations [7-10] or dissipative BBM equations [2,18]. This article is orga- 
nized as follows; in the next section we introduce the mathematical frame- 
work that we have chosen to study this dynamical system. In a third section 
we address the initial value problem for the evolution equation. In a fourth 
section we prove the existence of a smooth finite dimensional attractor. 

2. Mathematical framework 

2.1. Initial data. For the sake of convenience, we are interested in consid- 
ering periodic boundary conditions. We now consider functions for x G [0, 1] 
that are 1-periodic. We also assume that J f(x)dx = and / G L 2 (0, 1). 
Introducing w(t,x) = 1 + r)(t,x), we rather use the following system 

f ih ~ utxx ~ u xx + w x + u x u = f for (x, t) G (0, 1) x R + , , , 
\ w t + (wu) x -w xx =0 for (x, t) G (0, 1) x R+ . 1 ' 



The natural space for the velocity u(t, x) is then 

Hp er = {v£ H l per / J v{x)dx = 0}. (4) 

We then assume that no G Hp er . 
We now proceed to the assumptions on wq(x). 
The first physical assumption is that 

inf Wq(x) > 0. (5) 

This assumption ensures that the top of the fluid does not hit the bottom 
of the channel. The second assumption is that 

/ 1 ^o = l, (6) 
Jo 

that describes that w$ fluctuates around 1 the height at rest of the fluid, 
and that the fluids has constant volume. 
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The third assumption is related to the very definition of the entropy for 
convection equation, see [16]. 

Introducing Q(y) = ylny — y+1 that is convex and non negative, we assume 
that 

/ Q(w (x))dx < +00. (7) 
Jo 

Remark 2.1. Assume that f = here and that we are given regular enough 
solution (u, w) to (2) such that w > Vx, t . Multiplying (2) by (u, 1 + lnw) 
and summing the two resulting equations, we thus obtain 

d .t ,1 ,1„..2 



-\u\ 2 H1 + wlnw + \ux\l2 + / — = 0; 
2 JO Jo w 



(8) 



(n 2 + u 2 x )+ Q(w) 







rft L2 1 

t/ien, using' Jg 1 ttj = 1 

lim — 

i— >oo L2 

and i/ie /?md converges to the equilibrium (u,w) = (0, 1). 
2.2. Functional Analysis. We set 

K = {w(x) > 0; / w(x)dx = 1 and / Q(w(x))dx < 00}; (9) 

JO JO 

one may wonder which kind of topology we shall use in K. First of all, we 
observe that K is a convex set (this is obvious since the map y — ► Q(y) is 
a convex function). Furthermore, K is related to the following Orlicz space 
(see [1]). Introducing 

H(y) = (l+y)ln(l+y)-y, (10) 

that is a convex function, we observe that H and Q are related by the 
following inequalities. 

Lemma 2.2. 3Co > such that Vy > 

Q(y) - 1 < H(y) < C (Q(y) + 1). (11) 

Proof: On the one hand, if y G [0, 1], then Q(y) < 1 - y < 1 + #(y). If 
y > 1 then 

Q{y)=y\ny+l-y< (y + 1) ln(y + 1) - y + 1. 

To establish the reverse inequality, we observe that ylny ~oo (1+y) ln(y+l), 
then, for large y's , say y > R, H(y) < 2(nlny — y + 1). 
For y e [0,R], H(y) is bounded by H(R). Then the proof of the lemma is 
completed. □ 

Therefore w > belongs to K iff Jq 1 -u; = 1 and G L#, the Orlicz space, 
whose norm is defined by 

f 1 w 

\\w\\ Lh = inf{X >0, / H(-)(x)dx<l} (12) 
jo A 
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We now endow K with the topology of Lh , that is given by the distance 

d(wi,w 2 ) = \\wi - w 2 \\l h - (13) 

Remark 2.3. K is not a closed subset of Lh, but K = {w > 0; Jq w = 

1 and Jq 1 Q(w) < oo} is. Consider w n in K such that \\w n — w\\l h — ► 0. 
There exists \ n — > such that 

f 1 H( \ Wn ~ W \ )(x)dx< 2. (14) 
Jo \i 

We shall use 

H(y)> l -{^Tl-l) 2 . (15) 
Then, setting v n (x) = X~ 1 \w n (x) — w(x)\, 




= f(^+l-l) 2 + 2f(A^l-l) 
J o Jo 

< 4 / H(v n )+4([ (Vv~+T - 1) 2 )5 < 4(1 + ^2). 

JO JO 

XYien w n ^ w in L , and, up to a subsequence extraction w n — > u; a. e. □ 

3. The Initial Value Problem 
3.1. Main Theorem. 

Theorem 3.1. Consider the initial data (uq,wo) in H\ er x K. Then there 
exists a unique solution for J3J) 

(u(t),w(t)) G C(R+;H£ er ) x C(R + ;K) 

that satisfies moreover y/w G Lf oc (W + ; Hp er ) . 

Remark 3.2. In the theorem above, we would like to point out two important 
facts: 

* If inf wo (x) > 0, then, for any t,x > 0, w(t,x) > 0, that is physically 
relevant. 

* The dissipative Boussinesq system provides a smoothing effect in the w 
variable. In fact y/w belongs C(M+ — {0}; H%) 

Proof of Theorem 13.11 

Existence: we first regularize the initial data (u q ,Wq) to construct smooth 
solutions (u £ (t), w £ (t))'m C([0,T], Hp er x Hp er ) for instance. We then prove 
some a priori estimates and finally pass to limit. Since these methods are 
classical, we just indicate below how to derive the a priori estimates, refer- 
ring the reader to [3], for details. For the sake of simplicity, we drop the 
subscript e throughout the proof of the theorem. 
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First step: we first prove w(t, x) is positive . 
Consider a = inf wq > 0. Introduce 

J(t, x) = max(0, a. — w(t, x)). (16) 

Using the Kato's inequality ( see [12, 13] ) that reads (in the distribution's 
sense) 

w xx sgn(w) < (\w\) xx , (17) 

we thus obtain 

Jt + (uJ) x - J xx < 0. (18) 
Then, integrating in x, we have that 

J(t,x)dx<[ J(0,x)dx = 0, (19) 



and w(t,x) > a a.e. 

This result is related to parabolic Harnack's inequalities, see [6]. 

Second step: a priori estimate in Bl er x K. 
We begin with a technical lemma 

Lemma 3.3. Consider w > a smooth periodic function such that Jq 1 w{x)dx 
1 then 

W T. \ 2 



Q(w)dx< / ^V. (20) 
v Jo w J 

Proof: since — In is convex, Q(w) < w(w — 1) — w + 1. Then, using once 
more J w = 1, 



I Q{w) < f 1 w 2 -l. (21) 
Jo Jo 



On the other hand, for any ip smooth periodic function, then 

V\x)<2{j\ 2 )\j\ 2 f+j\ 2 . (22) 
We apply this to (p = \fw. Then 

jy< Mh ^j\)<{fjAy^ ,23, 

this concludes the proof of the lemma □. 

We now proceed to the a priori estimates. We multiply by (u, l+lnu>). 
We integrate in x over [0, 1] the resulting equations and then sum to obtain 



d_ 

dt 



1.. Il2 

w In w — w + 1 ) + 



+ J (uw) x In w 

+ / / ^ + 11^111,2= / M24) 

JO w Jo - j Q 



Integrating by part, we observe that 

(uw) x In w + / WiJi = 0. (25) 

o 
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On the other hand, using Young's and Poincare-Wirtinger inequalities, we 
obtain 



/ fu + \\u x \\ 2 L 2 
Jo 



— \\ u - 

1 

2 



z|lL 2 



2 



(*-\)\W\h 



2 

L 2 ' 



1 



47T 



2 

L 2 



47T 



We thus obtain, 
d - * 



dt 



Q(w) + -\\u\\ m 



+ ^\\u\\m 



THk < 1 

W ~ 47T — 2 



2 

L 2 ■ 



We now infer from (|21|) and (|23j) that 



W 



+ 



(26) 
(27) 
(28) 



We combine this inequality together with (f2~H) , we integrate with respect to 
t (thanks to the Gronwall lemma), and thus obtain 

Jl QM*)) + \\\u{t)\\ 2 m <\ + i^H/Hi 2 + e ~*(jf OW + ^IKHhO- 

(29) 

Since i7 is a convex function, for A > 1 (observe H(0) = 0), and due to 
Lemma 12.21 



H(w)dx < 



A 



1 + sup 

t>o Jo 



Q(w(t)) 



< 1, 



(30) 



for A large enough. Then w(t) remains bounded in the Orlicz space L#. 
We then have established an a priori estimate for (u, w) in L°°(i? + ; Hp er ) x 
L°°(R + ; L H ). 

Remark 3.4. Actually (27) implies that y/w is a.e. in t in H\. Since this 
Sobolev space is an algebra and since we can solve the evolution equation 
under consideration with initial data w in H\, this implies that for all t > 
w is in (smoothing effect). We precise this fact below. 

Third step: smoothing effect; ^Jw belongs to H 1 for t > 0. 
We set v = \/w that solves 



< 1 

v t -v xx h -u x v + uv x 

v 2 



0. 



(31) 



Multiply (USD) by —v xx and integrate. We then get 



1 d 2 

2Jt llVxllL2 + 



v 



xxWjJi 

1 

— 2 1 \^' xx 



h+c 



oo\\ u x\\L 2 



U X VV XX + 
2 



+ M 



UV T V 7 



oo\\ V x\\ L 2 



Then, since ||v||x, 2 = 1, 
d 



(32) 
(33) 
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Due to Gronwall lemma and since u is bounded in H 1 we then get the H 1 
bound on v = yjw. 

Remark 3.5. Actually for T < +00 y/w is in L 2 (0,T; H~) and then w in 
L\Q,T;Hl). 

Fourth step: uniqueness 

Consider two trajectories (u2>^2) and (ui,wx) that start from the same 
initial data. Due to the previous estimates both (112,^2) and (ui, w±) remain 
bounded in L 1 (0, T; ) . We set u = 112 — u\ , w = W2 — w\ that are solutions 
to 

ut ~ u txx - u xx + w x + \{u\ - u\) x = 0, , 34 s 
wt - w xx + (u 2 W2 - UxWx) x = 0. 

Then multiply these equations by (u, w) and integrate the resulting equation 
to obtain 

■^■(II^II^+II^IIl 2 ) ^ C(l+max(||u 2 ||Loo JI^IIloo, ||wi||l^)(||u||//i + ||w||l2). 

at 

The results follows promptly. □ 

4. The global attractor 
4.1. Existence of the global attractor. To begin with we state and prove 



Proposition 4.1. The semigroup S(t) defined on Hp er x K possesses 
absorbing set that is bounded in Hp er x H 1 



an 



Proof: The existence of a bounded absorbing set in Hp er x K comes from 
the estimate (|29|) of the previous section. Let to be the entrance time into 
this absorbing ball. Going back to ()33|) and applying the Uniform Gronwall 
Lemma (see Lemma III. 1.1 in [17]), we thus obtain that for t > 0, for some 
numerical constant c, 

t\\(V^)x(t + t )\\h <c(l + t)(l + ||/||l 2 ))exp(c + C ||/||| 2 ). (35) 

Therefore ^/w is bounded for large times into H 1 . Since H 1 is an algebra, 
then w is also bounded for large times in H . □ 

Theorem 4.2. The semigroup S(t) possesses a global attractor A in Hp~ er x 
Lh, that is a compact subset of H 2 x H 2 . 

Proof: we introduce the splitting (u,w) = (u 1 ,^) + (u 2 ,0), 
where u 1 satisfies 

u t ~ u txx ~ u xx + W x + U X U = / fQC\ 

u 1 (0) = 0, {6b) 



8 



M. ABOUNOUH, A. ATLAS, O. GOUBET 



and u 2 is solution to 

u t ~ u txx ~ u xx 

u 2 (0) = u . [6n 

We now define the families {Si(£)}t>o and {5*2(i)}t>o of maps in H 1 x Lh, 
where S\(t)(uo,wo) = (n 1 ,^) and Si(t)(uo,wo) = (u 2 ,0). 
First step : we prove that u 1 is bounded in H 2 . For this we multiply (|36j) 
—u xx and integrate between and 1 to obtain 



"XX 

1 d 



2dt llUxllHl 



p + IkLllia - - / fu l xx + J w x ii xx + J uu 
due to Young and Cauchy-Schwarz inequalities, then 

j t W u l\\m + \\<x\\h < c(\\f\\h + + IKII^IMIl-). ( 38 ) 

due to Proposition 14.11 l|29|). Gronwall and Poincare inequalities, we obtain 
that it remains in a bounded set of i? 2 for large times (t > to the entrance 
time into the absorbing ball). 
On the other hand, it is an exercice to prove that 

u 2 (t) — > strongly in H 1 when t — > oo. (39) 

Then Si(t)(uQ,wo) is bounded in H 2 x H 1 then compact in H 1 x Lh and 

S2(t)(uo,wo) — > in H 1 x Lh uniformly on bounded sets. 

Then from Theorem 1. 1.1 in [T] we have the existence of a global attractor 

A in H 1 x Lh, that is moreover a bounded set in H 2 x H 1 . 

We now prove that for a trajectory (u, w) in the global attractor, then w 

remains bounded in H 2 . For that purpose, multiply the second equation in 

© by W4 X and integrate by parts to obtain 



d 1 1 ll2 i, ,, 2 f 1 n f 1 1 f 1 2 

-dl^zll^ + ||uw|| L 2 = / u xx ww 3x + 2 u x w x w 3x -- u x w 2x 
at * jo jo ^ jo 



< c\\u\ 1^2 Hwll^i + -HiOotxII^. 

(40) 

Then the results follows promptly. It remains to prove that the global at- 
tractor, that is bounded in H 2 x H 2 , is in fact a compact subset of this 
space. This can be performed by the Energy Equation Method of [14] that 
is a suitable adaptation of the famous J. Ball argument. This is standard 
and will not be reproduced here; we refer the reader to [3] for details. □ 

4.2. Dimension of the attractor. In this section we are going to prove 
that the global attractor A has a finite dimension in £ = H 1 x {w G 
L 2 ;J w = 1}. £ is an affine space whose associated vector space is E = 
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H 1 x L 2 . To begin with, we need a result on the differentiability of the semi- 
group S(t) on the global attractor. Consider the non-autonomous linearized 
system 



v t - v txx ~ v xx + h x + (uv) x = 
h t + (uh + vw) x -h xx = 



(41) 



where (u(t),w(t)) = S(t)(uo,wo), (uq,wo) £ £, is a trajectory solution of 
and (wo,/io) £ E. Actually the linear mapping DS(t)(uo,Wo)i v o>ho) = 
(v(t),h(t)) is the uniform differential of S(t) as stated below 

Theorem 4.3. The non- autonomous PDE \4-l\j provides a well posed initial 
value problem in E. Moreover for T > 0, (vo, ho) £ E, (uq, wo) £ A, t < T 
there exists a constant C = C{T) such that 

\\S(t)(uo+vo,wo+ho)-S(t)(uo,wo)-DS(t)(uo,w )(v , h )\\ E < C(T)\\(v , h ) 

(42) 

where 1 < 5 < 2. 

Proof: to prove that the initial value problem is well-posed is standard 
and then omitted. Consider the solutions (u\(t), wi(t)) = S(t)(uo,wo), 

(u 2 (t),w 2 (t)) = S(t)(u +v ,w +h ) and (v(t),h(t)) = ^DS(t)(u ,w )^(v , h ). 

Then (p, q) = (u2,w 2 ) — (ui,w±) — (v,h) satisfies the system 



\ 6 E 



Vt ~ Ptxx ~ Pxx + Qx + {\v 2 +Vp + Uip+\p 2 ) x = 

Qt - q xx + (pq + ph + vq + hv + qu\ + w\p) x = 



(43) 



We shall use in the sequel that Jq 1 p = Jq q = and then ||p||#i and 
||Px||l 2 define equivalent norms. Multiply (|l3*|) by (p,q) and integrate to 
obtain (due to straightforward computations) 



ld_ 

2dt 



|2 , 
lL2 + 



2 

H 1 



+ IWxWh + \ \PxWh 



I ' qxP+l^ jip + v) 2 p x + J uipp x + J \pq + ph + vq + vh + wip + qui)q, 

< i\\v\\m + IWIl^OIMIhi I|9xIIl2 + WpxWl* 



+ 



(1 + \\ui\\ H i + \\v\\ H i + ||^l|| L 2 + \\h\\ L 2) 



1 1 2 . 1 1 1 1 2 
1 9a: 1 1 £2 + l|Px|li2 



+ 



L 2 \\qx\\ L 2 



We thus obtain, using the bounds on the attractor and the local in time 
bounds on (y, h) 



If 



\\{p,q)\\ 2 E < Ki(T)\\(p, q)\\ E + K 2 (T)\\(vq, ho)\\% + K 3 (T)\\(p, q)\\%. 



(44) 
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Consider a given interval of time [0,T]. Set e 2 = K 2 (T)\\(v , ho)\\% that 
is small. Then <f>{t) = exp(-iFCi(T))||(p, q)\\%(t) satisfies the ODE 

<j)<K<f) 2 + e 2 , (45) 
supplemented with 0(0) = 0. Then E(t) < 2e if e is small enough. □ 

We now give the main result of this section 

Theorem 4.4. The fractal and Hausdorf dimension in £ of the attractor A 
are finite . 

Proof: set £ = (u,w), j3 = (v, h). Now we study the operators DS(t)£o that 
contracts the m-dimensional volumes in E. Let /?q, ■■■jPq 1 in E. We study 
the following quantities 

G m = \\pHt)A...AP m \\l = det 1 < iJ < m (p\t),pi(t)) , (46) 

where f3 l {t) = (DS(t)S,o)Po- The Gram determinant G m represents the vol- 
ume of m-dimensional polyhedron defined by the vectors (3 (t), f3 m (t). We 
will show that for sufficiently large m this determinant decays exponentially 
as t — > oo. 

We consider (3(t) = (D5(£)£o)A) solution of (JUJ), we multiply by /? = (v, /i) 
and integrate to obtain 

Id f l f 1 f 1 f 1 

nirll^lll + \\ v x\\ 2 L 2 + \\h x \\ 2 L 2 = / uhh x + vwh x + uvv x + hv x . 
£ dt Jo Jo Jo Jo 

Recall that (u, w) is a trajectory that belongs to the global attractor. Intro- 
duce M = c(l + 1 1/| || 2 ) that is the H 1 bound for u in the attractor (see ((291 
)). We do not want to use estimates that involve H 2 norms of u as 
We bound the fourth term in the r.h.s of Q47|) by jH/ixll^ + IMI^- The 
third term can be bounded as follows 

, f 1 , 1„ 

| / UVV X I < -|M|#i|MIl 2 IMIl°° < 
JO ^ 

cIHIhiIHIx^HwxII^ 2 < \\\v x \\ 2 L 2 + cM^ 3 \\v\\ 2 L 2. 

We now proceed to the first term as follows 

-l i 
uhh x \ < -\\u\\ H i\\h\\ 2 L 4 < 

c||n||H 1 ||/ i || 3 All^llK 4 <^ll^ll! 2 + cM8/3 |l^ll^- 
For the second term, we have 



o 



BOUSSINESQ SYSTEM 



11 



f 1 1 

| / Vwh x \ < -7\\h x \\ 2 L 2 + |M|r2|M|!,cx> < 

Jo 4 

+ gll«*llia +c|Hlia|Hlia- 

To go further, we need a new estimate on w that reads 

Lemma 4.5. For any (u,w) in A, then \ \w(t)\\^2 < c(l + M 2 ). 

Proof: for a given trajectory in the attractor multiply the second equation 
by w and integrate to obtain 

r1 r 1 

IX iiiiO 11 11O I 11 11O 1 1 1 1 1 1 1 1 / 

-t\\ w \\l 2 + 2\\w x \\ L 2 = 2 uww x < \\w x \\ L 2 + \\u\\ La a\\w\\L°°', (48) 
at Jo 

here we have used that JqW = 1. We then infer from (|48|) that, using 
Poincare-Wirtinger inequality, 

/ (™-l) 2 = Hli*-l< IKIlia, ( 4 9) 
JO 

that 



d 1 

j t \\w\\ 2 L 2 + -\\w\\ 2 L 2<c(l + M 4 ). (50) 



Then the classical Gronwall lemma leads to the result. □ 
We then have 

1 d 

2dt' 

We introduce the Gram determinant 



+ ll« = c ((1 + M 8 )!!/^^) . (51) 



G m (*) = deh (A(^(t),^'(t))) 



where A (a, 6) = H a+6 Hg H a anc i that represents the m-dimensional vol- 
ume. Then we can proceed as in [7, 17] to establish that 

^ + mG m <c(l + M 8 )(y max min ^j^* 6 ' 1 ) G m . (52) 

Since the eigenvalues of the Laplace periodic operator are 47r 2 /c 2 each of 
multiplicity 2, then 

m \\R\\ 2 . m / 2 



V max min "" Tou^ - 2vr 2 V (2^)' 2 < (53) 

Therefore for m > c(l + M 8 ) the m-dimensional volume G m decays and the 
attractor has finite dimension. □ 
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